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Abstract 

We propose a stochastic particle model in (l+l)-dimensions, with one dimension corresponding to ra- 
pidity and the other one to the transverse size of a dipole in QCD, which mimics high-energy evolution 
and scattering in QCD in the presence of both saturation and particle-number fluctuations, and hence of 
Pomeron loops. The model evolves via non-linear particle splitting, with a non-local splitting rate which 
is constrained by boost-invariance and multiple scattering. The splitting rate saturates at high density, so 
like the gluon emission rate in the JIMWLK evolution. In the mean field approximation obtained by ig- 
noring fluctuations, the model exhibits the hallmarks of the BK equation, namely a BFKL-like evolution 
at low density, the formation of a traveling wave, and geometric scaling. In the full evolution including 
fluctuations, the geometric scaling is washed out at high energy and replaced by diffusive scaling. It is 
likely that the model belongs to the universality class of the reaction-diffusion process. The analysis of 
the model sheds new light on the Pomeron loops equations in QCD and their possible improvements. 
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1 Introduction 



Much of the recent progress towards understanding the dynamics of QCD at high energies 
comes from the observation [1,2] that the QCD evolution (at least, in its 'leading-logarithmic 
approximation' with respect to the energy logarithm In s) is a relatively simple classical stochas- 
tic process in the universality class of the 'reaction-diffusion process'. 

In its canonical formulation (see, e.g., Refs. [3,4]), the reaction-diffusion process A ^ AA 
involves particles of type A located at the sites of an infinite, one-dimensional, lattice, which 
are endowed with the following dynamics: a particle can locally split into two (A — > AA), 
two particles can recombine into one (AA — > A), and a particle can diffuse from one site to 
the adjacent sites. But the associated universality class covers a wide variety of stochastic pro- 
cesses, whose detailed microscopic dynamics can be different from those of the 'canonical' 
realization described above, but which share the same basic ingredients, leading to an effective 
one-dimensional dynamics with local growth and diffusion in the dilute regime, saturation of 
the occupation numbers at high density, and fluctuations in the particle number. These ingre- 
dients appear to be sufficient to ensure a universal behaviour for all such processes in specific 
limits (large time, weak noise), as empirically demonstrated by the experience with a large 
number of models [4], and conceptually understood from general arguments [5,6]. 

In particular, in the context of high-energy QCD, the 'particles' correspond to gluons (or 
dipoles), the 'time' axis corresponds to the rapidity y ~ In s, and the one-dimensional 'spatial' 
axis is the logarithm of the gluon transverse momentum (or the dipole transverse size). Further- 
more, the 'particle splitting' corresponds to the BFKL evolution [7,8], the 'recombination', to 
the non-linear effects responsible for gluon saturation [9-16], and the 'diffusion', to the non- 
locality of the BFKL kernel and the various gluon, or dipole, vertices in the transverse space. 
Of course, the actual QCD problem is much more complicated than just particle splitting and 
merging: it involves color degrees of freedom, whose role appears to be irreducible at high den- 
sity; also, gluon saturation (as described by the JIMWLK equation [14-16]) does not occur via 
'mechanical' recombination, but rather as a consequence of coherent, strong field, effects [17]; 
besides, the transverse plane has two dimensions and, moreover, the theoretical descriptions 
typically involve two types of transverse coordinates — the transverse size (or momentum) and 
the position in impact parameter space — , so the relevant transverse phase-space is truly four- 
dimensional. In spite of this, one can argue that, under suitable approximations, the dominant 
high-energy dynamics at fixed impact parameter is indeed one-dimensional, with the relevant 
dimension being the transverse size (or momentum), as alluded to above [1,2]. 

The universality of the reaction-diffusion process is especially interesting in this QCD con- 
text, since in spite of considerable work and progress over the recent years [1, 2, 12-27], a 
complete theory for the high-energy evolution and scattering in QCD is still lacking. More- 
over, even the approximate versions of this theory that are currently known, like the 'Pomeron 
loop' equations proposed in Refs. [2,20,21,24], or the 'self-dual' [23] effective Hamiltonian of 
Refs. [26, 27], appear to be too complicated to deal with in practice, and may also suffer from 
internal inconsistencies [28]. In view of this, it is particularly rewarding that some fundamental 
aspects of QCD at high energy, like the asymptotic behaviour of the scattering amplitudes (in 
particular, the breakdown of 'geometric' scaling [29-33] and the emergence of a new, 'diffu- 
sive', scaling at very high energies), could have been inferred from the correspondence with the 
reaction-diffusion process in statistical physics [1, 2, 6, 34-37]. 
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But this correspondence has its own, rather strong, limitations, which could be overcome 
only through deeper investigations of the actual, microscopic, dynamics. In view of the diffi- 
culty to perform such studies directly in QCD, one nowadays assists at a regain of interest in 
stochastic particle models [34,35,38^3], which aim at simulating the 'Pomeron loops dynam- 
ics' — by which we main the high-energy evolution and scattering in the presence of both 
saturation and particle-number fluctuations — in lower dimensions (zero or one transverse di- 
mensions). Such models are interesting not only because they allow for explicit solutions (even 
analytic ones, in the case of zero dimensions [40-43]), and thus permit us to familiarize our- 
selves with the physical consequences of the Pomeron loops, but also because they provide 
explicit evolution equations, and thus give us some insight into the structural aspects that one 
should expect in QCD. Besides, if sufficiently rich, a model may also allow us to study non- 
universal aspects, like the influence of the initial conditions at low energy, or the details of the 
transition from geometric to diffusive scaling. 

However, in order to meet with these desiderata a model should contain as much as possible 
of the actual QCD dynamics, to the extent that this can be made consistent with the limitations 
inherent in the model (like its low dimensionality and the lack of color degrees of freedom). 
For instance, the model should respect general symmetry properties, like the boost-invariance 
of the evolution equations, which in fact is one of the main physical arguments in favour of 
Pomeron loops [19,23,24,44]. (This condition is sometimes reformulated as the 'self-duality' 
of the evolution Hamiltonian [23,24,26].) Also, the model should incorporate some general 
properties of QCD at high energy, like the fact that a particle which scatters off a dense system 
undergoes multiple scattering. Moreover, in the limiting cases where the QCD equations are 
presently known, the model should reproduce suitable versions of these equations: it should 
be similar to the QCD dipole picture [8, 44] in the dilute regime, and it should reduce to a 
suitable version of the BK equation, or of the Balitsky-JIMWLK hierarch>Q [12-16], in the 
high-density regime where the particle-number fluctuations become unimportant. 

Together, the above conditions turn out to be very constraining, and it is not so easy to 
construct a model which fulfills all these constraints. In fact, most of the models proposed so 
far are directly inspired by the canonical reaction-diffusion process [34,35,38^1,43] (which 
is also closely related [43] to the old 'Reggeon field theory') and as such, not only they cannot 
be used to study deviations from universality, but also they have difficulties to accommodate 
multiple scattering: to be consistent with boost-invariance, such models must assume that a 
projectile particle can undergo only single scattering, however dense the target is. 

A noticeable exception in that sense is a zero-dimensional model originally proposed by 
Mueller and Salam [45], in an attempt to include saturation effects in the dipole picture, which 
has been recently reconsidered in Ref. [42] (see also Ref. [46]). In this model, boost-invariance 
and multiple scattering (in the eikonal approximation) are used as building principles, and it 
turns out that — due to simplifications associated with the lack of transverse dimensions — 
these two conditions are almost enough to uniquely fix the evolution law. The only additional 
assumption, which is fully consistent with the known situation in QCD, is that the evolution 
proceeds via particle splitting alone (no recombination), with one new particle being emitted 
per unit rapidity. Remarkably, the particle emission rate which emerges from these constraints is 



^ One should perhaps recall at this point that the BK equation and the Balitsky hierarchy are essentially 
equivalent at large-A^c> which is the limit of interest for the particle models. 
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non-linear in the number of preexisting particles and saturates at some maximal value when this 
number is large enough, as a consequence of multiple scattering. These features are appealing 
since rather similar to gluon saturation in the framework of the JIMWLK evolution. 

But the lack of transverse dimensions greatly restricts the spectrum of physical problems 
that can be studied in the framework of the model: the only non-trivial problem which can 
be addressed in zero dimensions [42,45] is that of the approach of the S'-matrix towards the 
black disk limit S* = 0. But the most interesting problems in the context of the high-energy 
evolution are those related to the coexistence of two regimes, one dense and one dilute (at low, 
and respectively, high transverse momenta), separated by the saturation momentum. In order 
to consider such problems one needs at least one transverse dimension — that expressing the 
gluon momentum, or the dipole size. It is our purpose in this paper to present an extension of 
the particle model in Refs. [42,45] to one transverse dimension, denoted as x, to be eventually 
identified — in the correspondence with QCD — to the logarithm of the inverse size of a dipole: 
X ln(rQ/r^), where r is the generic dipole size and tq is some arbitrary scale of reference 
(typically, the size of the dipole which started the evolution at F = 0). 

The model that we shall consider is constrained by the same physical assumptions as in 
Refs. [42,45] — boost-invariance, multiple scattering and evolution via particle splitting — , 
but in the presence of the transverse dimension these assumptions are not sufficient anymore 
to completely fix the dynamics. Yet, as we shall see, they strongly constrain this dynamics: the 
only non-trivial (class of) solution(s) to these constraints that we shall find is the one which 
appears as the natural, although not necessarily obvious, generalization of the original model 
in zero dimensions. Specifically, the elementary particle-particle scattering amplitude is now 
non-local in x, although very short-ranged, corresponding to the fact that QCD favors the 
interaction between dipoles with similar sizes. Accordingly, the rate for particle emission at a 
particular site x (i.e., for emitting dipoles of a given size r) is found to be proportional to the 
scattering amplitude between the new particle and the preexisting ones, and thus it depends upon 
the occupation numbers n{y) at all the other sites y. Once again, this emission rate saturates at 
its maximal value when the occupation numbers at the site of interest, or in its vicinity, become 
large enough. 

The one-dimensional model endowed with this dynamics develops a very interesting struc- 
ture, which is formally similar to that of the corresponding model in zero dimensions [42], but 
whose physical consequences are considerably richer, and also much closer to QCD. The evo- 
lution Hamiltonianffi which is self-dual (as it should, given the built-in boost-invariance), can 
be recognized as an extension of the JIMWLK Hamiltonian which allows for particle-number 
fluctuations. The equations for the scattering amplitudes generated by this Hamiltonian provide 
an intuitive generalization of the Balitsky-JIMWLK hierarchy in which the projectile and the 
target are symmetrically treated. The general structure of these equations is in fact similar to 
that of the 'Pomeron loop equations' proposed in QCD at large [2, 20, 21], in the sense of 
including three types of terms: linear terms responsible for the BFKL evolution, source terms 
describing particle-number fluctuations in the target (or, equivalently, saturation effects in the 
projectile), and non-linear terms corresponding to saturation effects in the target, which ensure 
the unitarization of the scattering amplitudes at high energy. 



^ In the context of a stochastic particle model, a Hamiltonian can be identified by rewriting the master 
equation in operator form; see Sect. |2]below for details. 
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At the same time, some interesting differences are found between the structure of the fluctua- 
tion terms in the toy-model equations and, respectively, the 'Pomeron loop' equations in QCD, 
which invites us to a more careful comparison between these two sets of equations. This study 
(to be detailed in the Appendix) shows that the observed differences reflect different ways to 
organize perturbation theory in the weak scattering regime, which are essentially equivalent, in 
the sense of providing the same dominant behaviour at high energy. This being said, the toy- 
model equations appear to be more symmetric in their treatment of the target versus projectile, 
and it is likely that a similar symmetry should hold in the complete equations of QCD, for which 
the toy model provides some inspiration. 

Furthermore, the presence of a transverse dynamics in the toy model allows for explicit 
comparisons with the corresponding dynamics in QCD and, more generally, for studies of the 
universality. It appears that the linearized version of the toy-model equations (as appropriate 
in the dilute regime) is very close to the BFKL dynamics, whereas the mean-field version of 
the same equations (in which the expectation values are assumed to factorize) is similarly close 
to the BK dynamics. Namely, in the dilute regime, the solution shows the characteristic BFKL 
pattern [7] (exponential increase in Y, exponential decrease in x with some 'anomalous dimen- 
sion', diffusion in x with diffusive radius oc VY), and even approaches 'color transparency' at 
very large values of x. After including non-linear terms in the mean-field approximation, one 
observes saturation, the formation of a traveling wave, and geometric scaling — all that being 
indeed very similar to QCD [30,31,33]. 

A stochastic particle model with the linear and mean-field behaviour aforementioned is ex- 
pected to belong to the universality class of the reaction-diffusion process, and this will be 
indeed confirmed by a numerical study of our model: the position of the traveling fronts (the 
toy-model analog of the logarithm of the saturation momentum) is thus found to be a random 
quantity with a dispersion increasing linearly with Y, and the shape of the average amplitude 
is seen to approach diffusive scaling at very large Y. Our present numerical analysis is only 
exploratory — it is merely intended to verify the universal features of the dynamics — , but we 
plan to perform a more detailed such a study in a further publication, and thus also study some 
non-universal aspects. 

To summarize, although conceptually rich and relatively close to QCD (at least, in the limit- 
ing cases in which the comparison is possible), the model that we shall introduce here remains 
relatively simple and allows for systematic studies, including numerical ones. This model could 
therefore serve as a playground for testing new ideas concerning either the formal structure of 
the high-energy evolution equations in QCD, or the effects of the evolution with Pomeron loops 
on various physical processes. 

The paper is organized as follows: In Sect.[2]we shall describe the construction of the model 
from the underlying physical assumptions, present the master equation for the probabilities, 
and discuss some of its general properties. In Sect. [3l we shall use the master equation to de- 
duce evolution equations for physical observables, like the A;-body particle densities and the 
scattering amplitudes. Then, in Sect. IH we shall present an analytic study of these equations, 
which will successively focus on the linear approximation, the mean-field approximation, and 
the qualitative aspects of the full evolution with fluctuations. These analytic considerations will 
be then illustrated by the numerical results in Sect. [51 which will also substantiate our expec- 
tation that the model should belong to the universality class of the reaction-diffusion process. 
Finally, Sect. [6] contains our conclusion and a brief account of the analysis in the Appendix. 
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2 Construction of the model 



In this section, we shall formulate our stochastic particle model for high-energy evolution 
and scattering in QCD. As explained in the Introduction, this is a (1 + l)-dimensional model, 
where one of the dimensions refers to the rapidity Y (the logarithm of the energy), which plays 
the role of a 'time' for the high-energy evolution, while the other one is a spatial dimension 
— the position of the particle along an infinite one-dimensional axis. For the analogy with 
QCD, one should keep in mind that this spatial dimension corresponds to a size in QCD, more 
precisely to the logarithm of the inverse of the transverse size of a dipole. For the sake of 
the presentation, we shall first describe the discretized version of the model, where the spatial 
dimension is assimilated to a one-dimensional lattice, with lattice points labeled by i. This is, 
of course, also the version of the model which is best suited for the numerical implementation 
to be eventually discussed. However, as we shall later demonstrate, the continuum limit of the 
model is well-defined and rather straightforward to obtain. 

A system of particles is then defined by specifying the number of particles rii (or 'occupation 
numbers') at all the sites of the lattice (some of these sites could be empty). With increasing 
rapidity, the composition of the system can randomly vary, via particle splitting according to 
a law to be shortly specified. Hence, the particle distribution is stochastic, and the state of the 
system at rapidity Y is described by the probability P{{n}, Y) to find a given configuration 
{n} = {..., rij, nj+i, . . . } of particles. 

In order to mimic a scattering problem, we shall consider two types of particles, say particles 
of type 'right' (R) and particles of type 'left' (L), which coexist at the lattice sites and interact 
with each other via generally non-local interactions (see below for details). The particle labels 
'i?' and 'L' refer, of course, to the fact that in an actual scattering problem, the two colliding 
systems propagate in opposite directions along the collision axio. 

In what follows, we shall assume that the average S-matrix element for the elastic scattering 
can be given the factorized form 

{S)y= Y1 PR{{n},Y-Yo)PL{{m},Yo)S{{n},{m}), (2.1) 

{n},{m} 

where Y is the total rapidity separation between the two systems and is divided between the 
'target' (the right mover), which has a rapidity Y — Yo'in the chosen frame, and the 'projectile' 
(the left-mover), which has a rapidity — Yq (with Yq > 0). Furthermore, P^^{{n}, Y — Y^ and 
Pl({'^}, ^o) are probability distributions describing the two systems at the time of scattering, 
and S(\n\, \j^Y) the ^-matrix for the elastic scattering between two given configurations 
{n\ and {m} of the target and, respectively, the projectile. Finally, the sum in Eq. (12.11 ) runs 
over all the possible configurations in both systems. 

A factorization similar to Eq. (12.11 ) is known to hold in QCD for onium-onium scattering 
within the dipole picture [44,47]. More generally, such a factorization emerges from a more 
fundamental, quantum, description in terms of light-cone wavefunctions whenever the single- 
particle states are eigenstates of the S-matrix operator (which is our implicit assumption here). 



° One should not confound this collision axis with the spatial dimension in the toy model, which rather 
corresponds to transverse degrees of freedom in QCD, as explained before. 
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On physical grounds, the average S'-matrix must be independent of the rapidity divider Yq, 
i.e., upon the choice of a Lorentz frame, which implies (with short-handed notations): 



= ^ = E [MY - Y,) + — PdYo) S{{n}, {m}) . (2.2) 

This condition, which relates properties of the scattering to those of the evolution with Y, rep- 
resents a rather strong constraint on the latter, and can be even used to fix the structure of the 
evolution law under some additional assumptions. Let us enumerate here the main assumptions 
that we shall rely on in that respect: 

( i ) Multiple scattering in the eikonal approximation. Let Cij = 1 — Tij denote the S'-matrix 
describing the scattering between two elementary particles with positions i and j ; here, 
Tij is the corresponding T-matrix and is assumed to be real, as expected for the dominant 
behaviour at high energy. Each particle from the target can scatter off all the particles in 
the projectile, and different particles from the target scatter independently off each other. 
Then, the S-matrix for given configurations of the target and the projectile reads: 

Si{n},{m}) = n^r^- (2.3) 

id 

(ii) One particle emission per evolution step. We shall assume that, when increasing rapidity 
in one step (Y Y + dY), the evolution consists in the emission of a single additional 
particle, with a probability which generally depends upon all the occupation numbers in 
the system. That is, the additional particle is emitted coherently from all the previous ones. 

Both these assumptions are natural from the perspective of perturbative QCD at high en- 
ergy: The eikonal approximation is the standard method to resum multiple scattering in QCD at 
high-energy, for both gluons (see, e.g., Refs. [12, 14, 15]) and dipoles [8, 13,44]. Furthermore, 
the high-energy QCD evolution in the leading logarithmic approximation consists indeed in 
the emission of a single s-channel gluon at each step of the evolution. Note that, unlike other 
models in the literature [39-41,43], our model does not include a mechanism for particle recom- 
bination, in agreement with the fact that there is no (s-channel) gluon, or dipole, recombination 
in QCD in the leading-logarithmic approximation. Rather, as we shall see, particle saturation 
occurs via high-density effects in the emission rate, similar to the gluon saturation in the frame- 
work of the JIMWLK evolution [17]. 

In zero (transverse) dimensions, that is, for a model living at a single lattice site, the as- 
sumptions (i-ii) above are sufficient to completely fix — via the condition (12.21) for boost 
invariance — the structure of the evolution equation for the probabilities P [42,45]. But with 
one spatial dimension at our hand, there is still much freedom left, because one can make dif- 
ferent choices for the non-locality of the splitting process. As we shall see in what follows, a 
relatively simple model which is conceptually interesting (in the sense of bearing many similar- 
ities to QCD) can be obtained after also imposing the following, third, constraint: 

( i i i ) After one step in the evolution, the configuration of the system changes only by the addition 
of one new particle at some arbitrary site. That is, the configuration remains the same as 
prior to the evolution except for the presence of the additional particle. This is clearly a 
stronger constraint than the above assumption (ii) which would allow for any final state 
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with one additional particle as compared to the original state, irrespective of the spatial 
distribution of the particles in the final state. 

Under this additional assumption, the evolution equation for P({n}, Y) (the 'master equation') 
takes the following generic form: 



dP{{n},Y) 
dV 



J2[M...,n,-l,...)P{...,n,-l,. ..,¥)- M{n}) P({n}, Y)] , (2.4) 



where the quantity fi{{n}) has the meaning of a 'deposit' rate: this is the probability per unit 
rapidity to find an extra particle at site i after one step of the evolution, starting from an original 
configuration {n}. The positive (gain) term in the r.h.s. Eq. (12.41 ) shows that in order to arrive 
at a final configuration {n} after one step in this evolution, one needs to start in a configura- 
tion containing one particle less on the site i, where i is arbitrary. The negative (loss) term is 
necessary for the conservation of the probability. 

The evolution law in Eq. (12.41 ) deserves some physical discussion. But let us first complete 
the construction of the model, by deducing an explicit form for the deposit rate fi{{n}). As 
anticipated, this is obtained by exploiting the condition of boost-invariance, that is, by inserting 
the generic master equation (|2.4I) into Eq. (|2.2I) . to deduce the following constraint : 

E if^iM) - /^( W) t^i{n})] = 0, (2.5) 

i 

where we have defined 

i^(w) = i-n^.7- (2-6) 

The above notation is not accidental, since ti({n}) is formally equal to the amplitude for the 
scattering of a projectile particle at site i off a target with a given configuration {n}. The most 
general solution to Eq. (12.51) that we have found reads 

= with Cij = Cji, (2.7) 

j 

and where Cij is independent of the particle occupation numbers. It is convenient to choose 
Cij oc Sij, since this leads to a relatively simple and intuitive model which still has a rather rich 
structure, as we shall later see. Namely, we shall fix our model by choosing 

A 

Cij = — Sij , (2.8) 

where r = th and A stands for the "lattice spacing" along the discretized spatial axis. The 
fact that Cij is proportional to A is natural, since the deposit rate density fi/ A will be well- 
defined in the continuum limit A — > 0. The fact that we choose Cij proportional to 1/r is again 
natural, since in the low density limit fi should be of order 0{1) when measured in units of r. 
Indeed, when the n^'s are relatively small, such that UjTij -C 1, the scattering amplitude tj in 
Eq. (12.61) reduces to U ^ J2j and therefore the particular parametric dependence of Cij on 
r exhibited above is clearly necessary for the condition to be fulfilled. 
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To conclude, as a result of the various assumptions presented above, we finally arrive at the 
model defined by the master equation (12.41 ) with the following expression for the deposit rate : 



/.(M) l-n.crg U{{n}) 

A T T 



(2.9) 



Up to an overall normalization, this is simply the scattering amplitude for the scattering between 
the particle newly produced at site i and the preexisting particles in the system, from which the 
new particle has been emitted. 

To better appreciate the physical content of the above formula, we need to know a little bit 
more about the elementary scattering amplitude Tij. Later on, we shall present an explicit model 
for this quantity, dictated by the analogy with QCD (see Sect. |4]for details). Here, it suffices 
to mention that Tij depends only upon the separation \i — j\ and is rapidly decreasing when 
increasing this separation. In particular, the quantity r = th is independent of i (a property that 
has been already used in writing Eq. (12.81) ) and will be taken to be small: r ^ 1 (the analogous 
quantity in QCD is of C(a^)). One therefore has <^ 1 for any pair (ij). 

Armed with this knowledge, we now return to Eq. (|2.9I) and discuss some limiting cases: 

(i) When the occupation numbers are relatively small, such that uj <^ l/rij, one can ex- 
pand a^j 1 — TijTij, and then the deposit rate density for the given site i becomes simply 
proportional to the particle occupation numbers at all sites: 

i^l^Jl ^Y.^^-j when rij < l/nj for any j. (2.10) 

This describes a situation where the extra particle at i is incoherently emitted by any of the 
preexisting particles in the system. When the approximation (12.101) holds for any site i, we shall 
speak about the dilute, or linear, regime. Since, as mentioned before, Tij is a decreasing function 
of the separation \i — j\, the dilute regime is realized when all the occupation numbers obey 
the condition rij <^\/ r. The dipole picture in QCD [8] is the right term of comparison for our 
toy model in this dilute regime; and, indeed, if one computes the analogous deposit rate in the 
dipole picture, one finds that this is linear in the dipole density [44]. 

(ii) When at least one of the occupation numbers rij becomes so large that a^^ <^ 1 (this 
requires nj ^ 1/ Tij), the deposit rate at i saturates at its maximal value of (9(l/r) : 

•^'^^ ^ — when rij ^ l/rjj for some j. (2.11) 

In particular, a sufficient condition to have saturation at a given site i is that the respective 
occupation number rii obeys rij ^ 1/r. This is again similar to the situation in QCD, where the 
gluon emission rate in the JIMWLK equation (the analog of the deposit rate in the toy model) 
saturates when the gluon occupation numbers are sufficiently high. 

One may also notice that in the limit where the scattering is truly local, tij = Tdij, Eq. (12.91) 
yields /j/ A = ( 1 — cr"' ) / r, with a = 1 — r . In this limit, the different lattice sites decouple from 
each other, and each of them evolves according to the zero-dimensional model of Refs. [42,45]. 
This limit is uninteresting for the present purposes, and shall not be considered anymore. 

Let us now return, as promised, to a physical discussion of the evolution described by 
Eq. (12.41 ). To that aim, it is useful to imagine a specific microscopic mechanism leading to 
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this evolution. It is natural to assume that the process leading to the production of the extra par- 
ticle in one evolution step is the splitting of one particle into two new ones. For this mechanism 
to be consistent with Eq. (12.41 ). one additional assumption is however needed: at least one of the 
two daughter particles produced after such a splitting must remain on the same lattice site as 
its parent particle. Introducing the emission rate gj^ji({n}) for a particle at j to split into two 
particles at j and i, respectively, with arbitrary i, we are lead to the following master equation 



dP{{n},Y) 



^ [gj^.l . . , - 1, . . . ) P(. . . , - 1, . . . , F) - g,^,^{{n}) P{{n}, Y) 



dV 

(2.12) 

which leads to Eq. (12.41 ) after identifying fi{{n}) = J2j The process in Eq. (12.121) 

is quite peculiar in the sense that the particle which splits does not disappear in the final state, 
rather it gets replaced by a new particle at the same site[^. In the analogy with QCD — in which, 
we recall, particle positions correspond to dipole sizes — , this would correspond to a process in 
which a dipole of size j splits into two dipoles, one with the same size j and the other one with 
some arbitrary size i. This is quite different from the (BFKL) picture of dipole splitting in QCD 
at large Nc [8], where the daughter dipoles have generic sizes, which typically are comparable to 
the size of the parent dipole. But in spite of this explicit dissimilarity at the level of the physical 
interpretation, the toy model evolution generated by Eqs. (12.41) and (12.91 ) turns out to have many 
common features with the high-energy evolution in QCD, as it should become clear from the 
analysis in the subsequent sections. 

Last, but not least, we notice that one can put the model in the Hamiltonian form 



dPi{n},Y) _ A 
dV r V 



exp 



d 



dm 



1 - exp ^ ^ Tij In cTij^ P{{n}, Y) 



= HP{{n},Y), (2.13) 

where the differential operator exp (—d/dui) is the 'translation' operator which reduces the 
occupation number at site i by one: exp (—d/drii) F{. . . , nj, . . . ) = F{. . . , nj — 1, . . . ) for a 
generic function F{{n}). This rewriting makes it explicit that the Hamiltonian underlying the 
evolution of the model is 'self-dual', that is, it is invariant under the self-duality transformation 
[23] which in the present context consists in exchanging 

d 

— < > -Y^rij In ai,, (2.14) 

* j 

and then reversing the order of the operators. This self-duality is the expression of the constraint 
of boost-invariance on the structure of the evolution Hamiltonian [23,24]. Thus, this condition 
is automatically fulfilled by our model, where boost-invariance is built-in. Also, the presence 
of two types of exponentials — one involving the particle occupation number and the other one, 
the derivative with respect to it — in Eq. (12.131) is reminiscent of the two types of Wilson lines 
which appear in the QCD Hamiltonian proposed in Refs. [26,27]. In that case, Y.j %| Incjjl 



^ Incidentally, this discussion also shows that a natural generalization of the present model would be to 
allow for the more general splitting i —>■ jk with generic i, j, and k, and with a rate gi^jk{{i^})- 
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is replaced by the color field produced by the s-channel gluons, and d/dui by the (functional) 
derivative with respect to the color charge density of these gluons. 

It is furthermore interesting to consider the approximate version of the master equation which 
is formally obtained by expanding the operator exp {—d/drii) to linear order in the derivative: 

^^^fe^ - -T.-^;W{n})P{{n}.Y)\ (2.15) 

Less formally, this is simply the limit in which the difference between the two terms in the r.h.s. 
of Eq. (12.41 ) is assimilated to a derivative with respect to rii, which is reasonable when rii ^ 1, 
for any i. This is the limit of large occupation numbers, in which the fluctuations associated 
with the discreteness of the particle number become negligible. 

In this limit, the self-duality of the evolution is of course lost, since the system property 
of being dense depends upon the frame. Moreover, it should be intuitively clear that, whatever 
frame we choose, the system cannot be dense at all the sites. Indeed, the spatial axis being 
infinite, there will always be an infinite number of sites where the occupation numbers are either 
zero, or small, of 0{\), and this for any Y . Hence, Eq. (|2.15l) is merely a formal approximation, 
which for a given Y has, at most, a limited range of validity in x — it refers to the evolution of 
the bulk of the particles in a region with high occupancy. 

This discussion is reminiscent of that of the JIMWLK equation in QCD, which has been 
established [14, 15] for systems with a relatively high gluon density and which ignores the fluc- 
tuations in the gluon number [2] . And, indeed, one can recognize the above equation (|2.15l) as 
the toy-model version of the JIMWLK equation. This analogy has been already discussed in the 
context of the zero-dimensional model in Ref. [42], to which we refer for more details. To sum- 
marize this section, the toy-model that we have here introduced reduces to (one-dimensional 
versions of the) dipole picture in the dilute regime and, respectively, to the JIMWLK picture 
in the regime of high occupancy, and in general it provides a self-dual 'interpolation' between 
these two pictures, as necessary for boost-invariance. 



3 Evolution equations for the observables 



Given the master equation (12.41) along with the deposit rate, Eq. (12.91) . one can obtain the 
evolution equation for any 'observable' O, by which we shall mean a quantity whose event- 
by-event value 0{{n]) depends upon the configuration {n} of the particles in the system. Its 
average value at rapidity Y , which is a measurable quantity, is given by 



{0)y = Y.P{{n},Y)0{{n}). 

{n} 



(3.1) 



Differentiating with respect to Y and performing a shift rii — l — in the contributions arising 
from the gain (plus sign) terms in Eq. (12.41) one easily finds 



dip), 
dY 



^(/,(M) [o(...,n, + i,...)-o(M) 



Y 



(3.2) 



11 



In what follows, we shall specialize the generic evolution equation (13.21 ) to two interesting types 
of observables, namely, the A;-body particle occupation numbers and the scattering amplitudes 
for projectiles made with a given number of particles. In the process, we shall also demon- 
strate that the continuum limit A — > of the ensuing equations is well defined, and we shall 
furthermore compare these equations to the corresponding ones in QCD. 



3. 1 Particle number densities 



The simplest observable that one can think of is (nj)y — the average occupation number at 
the site i. From Eq. (13.21 ). it is straightforward to obtain 

(3.3) 

with an obvious physical interpretation: the rate for the change in the average particle number at 
a given site is equal to the average value of the deposit rate at that particular site. Two limiting 
cases of Eq. (13.31) are particularly interesting, corresponding to the two limits of the deposit rate 
fi discussed in relation with Eq. (12.101) and, respectively, Eq. (12.1 II) : 

(i) If the system is relatively dilute 'around the site i' — meaning that, for the typical 
configurations, one has ^ 1/r and also rij <^ 1/r for all the sites j which are not too far 
away from i — , then we can linearize (/i)y with respect to the (average) occupation numbers 
and thus obtain: 

^-7E-.- (3-4) 

In this regime, the particle occupation at site i and at the neighboring sites undergoes a linear 
evolution, leading to a rapid (exponential in Y) rise in the respective average occupation num- 
bers. This rapid rise, which will be explicitly demonstrated in Sect.lH is the toy-model analog 
of the BFKL growth in the gluon (or dipole) distribution in QCD in the dilute regime [7]. 

(ii) If, on the other hand, the system is 'dense at i\ meaning that the average deposit rate 
(/j)y saturates at its maximum value, cf. Eq. (12.1 II ). then the growth of the (average) occupation 
number at i is considerably slowed down: 

^ « - ■ (3.5) 

oY T 

This equation shows that, once the particle occupation number at i becomes of 0{1/t), its 
subsequent growth with Y is only very slow, linear in Y. This is the mechanism of saturation 
within the toy model, and is similar to gluon saturation in QCD, where the occupation numbers 
for the saturated modes grow also linearly with Y [17,48,49]. 

Returning to Eq. (13.31) . one should also notice that this is not a closed equation for (nj)y, 
but rather the beginning of a hierarchy : if one expands a^j = exp{nj In aij} in powers of rij 
inside the brackets in the r.h.s., one generates /c-body correlations of the occupation numbers 
with any k > 1. This invites us to consider the evolution equations obeyed by the general, 
'normal-ordered', /c-body occupation numbers, defined as 
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nf2..ik = ^iMi2 - ^hi'z) ■ ■ ■ i^i, - Si,i^ ^h^iij- (3-6) 

Here, the 'normal ordering' refers to the subtraction of 5-functions, as explicit in the equa- 
tion above, which are needed to guarantee that, in constructing n^'^\ one counts only sets of k 
particles which are all different from each other. The corresponding evolution equation can be 
readily obtained as 

^%^=i:(/.,(W)"M„.j^, (3.7) 

where the notation ^ means that the ij variable is to be omitted. 

Before proceeding with the hierarchy for the scattering amplitudes, let us first briefly discuss 
the continuum limit of our model, which is straightforward. Clearly, this limit amounts to re- 
placing i Xi = iA and then letting A ^ 0, so that Xj x, with x our continuous spatial 
variable. Then, a sum over i gets converted into an integration over x: 

^Y.F^^ J dxF(x), (3.8) 

i 

while a sum over all the possible configurations of the onium wavefunction becomes a path 
integral, that is 

T.^ f [Dn{x)]. (3.9) 

{n} 

The /c-body (normal-ordered) particle densities are obtained from the corresponding occupa- 
tion numbers after dividing by A''' and then letting A ^ 0. For example, the 2-body, or pair, 
density is given by 

^ a" Z ) " ^'^^^ ^^^"^ - -^(a; - y) n{x), (3.10) 
while the deposit rate density can be immediately obtained from Eq. (12.91) . and reads 

X _ 1- MM) _ 1 -exp [/ dzn{z) \n(T{x\z)] _ t{x) 
f[x)-nm ^ - ^ _ ^ . p.iij 

It is now straightforward to reexpress all the equations presented earlier in continuum notations. 
For instance, the 1-body and 2-body particle densities obey the following equations, as obtained 
by taking the continuum limit in Eqs. (|3.3|) and (|3.7I) (with k = 2), respectively: 



d {n{x))y 

-\J K^J/Y^ 



{f{x))y, (3.12) 



d (n^'^\x, y) 



— ^ = {f{x)n{y) + f{y)n{x))y . (3.13) 
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It is important to notice that Eq. (13.131) has no singularities in the equal-point limit y = x 
thanks to the subtraction of the 5-function in the definition (13.101 ) of the normal-ordered 2- 
body density. A similar property holds for the A;-body densities with k > 2. 



3.2 Scattering amplitudes: Beyond the Balitsky-JIMWLK equations 

We now move to the scattering problem and start by assuming that the projectile consists in 
a single particle located at i. In QCD, this would correspond to a projectile made with a single 
dipole of size oc exp(— Xj/2). This implicitly means that we work in a frame in which almost 
all of the rapidity Y is carried by the target, whose wavefunction evolves according to Eqs. (12.41 ) 
and (12.91 ). whereas the rapidity of the projectile is so low that its evolution can be neglected. 
The ^-matrix describing this scattering is obtained by replacing mj 5jk in Eq. (12.31) . which 
then yields Sk{{n}) = 11/ cT/fc • By making use of Eq. (13.21) . we deduce the following evolution 
equation 

^ = (3.14) 

j 

which is not a closed equation — the average S-matrix element for the one-particle projectile 
being related to that for a projectile made with two particles — , but only the first equation in 
an infinite hierarchy. The general equation in this hierarchy can be obtained by studying the 
scattering of a projectile made with m particles, at given position^i] ii,i2, im- For a given 
target configuration {n}, the corresponding ^-matrix reads 

Sn■■■s^^ = l[^l■■■]l^^^ = ]l (^n. • • • ^^™.)"^ (3-15) 

jl jm j 

and then Eq. (13.21) implies 

d {si, . . . Sijy _ ^ A(l - a,,j . . . (Ji^j^ 



dY 



E ' "--^-----^^ . . . S^Jy . (3.16) 



Although obtained here by following the evolution of the target, these equations can be easily 
reinterpreted as describing evolution in the projectile. To that aim, it is important to notice that 
the kernel in the above equation, i.e. the quantity in front of the expectation value in the r.h.s., 
is precisely /j({m}), that is, the deposit rate at site j evaluated for the given configuration 
{m} of the projectile. Thus, Eq. (13.161) can be interpreted as follows: a splitting takes place in 
the projectile, leading to a system with m + 1 particles (with j being the position of the extra 
particle). Subsequently this new system scatters off the target giving rise to the first term in the 
above equation, while the second (negative sign) term there corresponds to the possibility that 
no extra particle was created. 

The physical content of this hierarchy will be further discussed on its continuum version, 
which reads (cf. Eq. (12.91) ) 

^ Of course, some of these m particles can have identical positions. 
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I ■ ■ ■ ^Xm) 

dY 



I mm}) (5: 




(3.17) 
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in simplified notations where the coordinates are shown as lower indices and the F-dependence 
of the expectation values is not indicated. The general structure of these equations is interesting, 
as it may shed light on the corresponding structure in QCD. It is especially instructive to com- 
pare these equations to the Balitsky-JIMWLK equations [12-16] and also to the more complete 
'Pomeron loop' equations recently proposed in QCD at large [2,20,21]. 

Recall at this point that the Balitsky hierarchy in QCD has been obtained by performing 
different evolutions in the projectile and the target: the projectile has been assumed to be dilute 
and evolve according to the dipole picturelZl whereas the target was taken to be dense and 
evolve according to JIMWLK equation. Under these assumptions, the Balitsky equations are 
formally boost invariant, but the whole scheme is clearly incomplete and unsatisfactory, since 
an evolving projectile eventually becomes dense and, vice versa, even a target which looks 
dense on some resolution scale (at relatively low transverse momenta) has necessarily a dilute 
tail at high transverse momenta, and this tail is the driving force for its evolution [1]. 

This problem is overcome, by construction, within the toy model, where the target and the 
projectile are symmetrically treated, and it is interesting to see how this is reflected in the struc- 
ture of the evolution equations. The first equation in the hierarchy (13.171) . that is. 



is formally similar to the corresponding Balitsky equation, from which it differs only by the 
replacement of the 'dipole kernel' (the rate for dipole splitting in QCD) by the reduced scat- 
tering amplitude t^z/t, which plays the role of the elementary splitting rate (corresponding to 
a parent particle in isolation) within the toy model. Note that the splitting is more constrained 
in the toy model than in QCD, for the reasons explained in Sect. |2] one of daughter 'dipoles' 
X and z which appear in the r.h.s. of Eq. (|3.18l) is bound to have the same 'size' x as its parent 
'dipole' in the l.h.s. 

It seems therefore natural to identify the elementary splitting x ^ xz with rate t^z/t in the 
toy model to the dipole splitting in QCD. With this identification, the differences between the 
toy-model hierarchy (13.171) and the Balitsky equations become visible in the higher equations 
in the hierarchy, starting with the second one: 



The respective Balitsky equation would involve only the two positive terms, r^^ and Tzy, in the 
numerator in the kernel, which describe the independent splittings of the particle at x and at 
y, respectively. The additional, negative, term —r^zTzy corresponds to 'saturation effects' in the 
evolution of the projectile, namely to the fact that the two splittings are not truly independent. 



^ We here restrict ourselves to the large- A'^c version of the Balitsky hierarchy, as appropriate for com- 
parison with the toy model. 




(3.18) 



z 




(3.19) 



z 
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This term is formally suppressed by a power of r with respect to the previous ones, but as we 
shall shortly argue it cannot be neglected, since it plays an essential role in the evolution. 

This discussion can be easily generalized to the higher equations in the hierarchy: the toy- 
model analog of the Balitsky equations can be obtained from the general equations (13.171) by 
replacing the kernel fz{{'m}) there by its linearized version (cf. Eqs. (12.101) and (13.111) ) 

mm}) — E , (3.20) 

i=i 

which describes the independent splitting of any of the m particles in the projectile. 

Clearly, this is similar to QCD, since it is tantamount to saying that the projectile evolves ac- 
cording to the dilute approximation (12.101) to the master equation (|2.4I) . Also like in QCD, it 
can be easily checked that the same set of equations would be obtained by evolving the target 
according to the toy-model version of the 'JIMWLK' equation, Eq. (12.151) . 

Thus, clearly, the toy-model hierarchy (13.171) goes beyond the Balitsky- JIMWLK equations 
by including saturation effects in the evolution of the projectile, or, equivalently, particle- 
number fluctuations in that of the target. It is therefore interesting to better understand the 
physical role of the additional term, so like the negative term in the kernel in Eq. (13.191 ). To 
that aim, and also to facilitate the comparison with the 'Pomeron loop' equations of QCD, it 
is convenient to rewrite the above equations in terms of the T-matrix elements (or 'scattering 
amplitudes') = 1 — s^- We display here the first two equations in the ensuing hierarchy (the 
third such equation will be shown in the Appendix): 



^^^^^ ^ ""^'^ (t. - M.) , (3.21) 



dY J T 

z 

and, respectively. 



dY 



— {{tz - U,)ty) + ^ ((t, - tyt,%) + - t,)(l - ty)) 

T T T 



.(3.22) 



Notice that t^ty = (1 — Sa;)(l — Sy) is the amplitude for the simultaneous scattering of two 
particles in the projectile. At this level, it becomes straightforward to recognize the toy-model 
analogs of the BFKL and, respectively, BK equation: the former is obtained by neglecting the 
non-linear term {t^tz) in the r.h.s. of Eq. (13.211) . while the second corresponds to treating this 
term in a mean-field approximation (MFA) which assumes factorization: {txtz) ~ {tx) {tz)- 
However, this factorization is inconsistent with the higher equations in the hierarchy, starting 
with Eq. (13.221) . This breakdown of the MFA is associated with the presence of 'fluctuation 
terms', like the last term oc r^zTzy in the r.h.s. of Eq. (|3.22l) . which are precisely the additional 
terms with respect to the (toy-model) Balitsky equations. Such terms are important, since they 
generate the correlations in the dilute regime, as we explain now: 

Consider indeed the regime in which the target is so dilute that the average scattering am- 
plitude is very wealj^: (t) < r. This corresponds to the dilute tail of the target distribution 
at large values of x, where the (average) occupation number is {n) < 1. (Indeed, in the dilute 



The relevance of this regime for the high-energy evolution will be explained in Sects. |4]and|5] 
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regime, one can write (t^.) ^ t^z {n^) ; see, e.g., Eq. (12.61) .) In this regime, the last term in 
the r.h.s. of Eq. (13.221) can be approximated as (T^zTyz/r) {tz), which for (t) < r is comparable 
to, or even larger than, the BFKL-like terms in that equation, so like (t^z/t) {tzty). This shows 
that, in this dilute regime, the two-particle amplitude (tt) is still very small but it gets built from 
the one-particle amplitude (t), via the last term in Eq. (13.221) . To better appreciate the physical 
interpretation of this term, notice that in the dilute regime one can successively write 



d {t^ty) 
dY 



fluct 



{tz) ^ / -^^^ I Tzy, (n^) ~ / T^rzTyz Qy , (3.23) 



where we have also used the appropriate limit of Eq. (13.121) . According to this equation, the 
change in the two-particle amplitude {txty) in one step of the evolution can be interpreted as 
the following evolution in the target: first, a new particle is created at z, via the splitting process 
w — i> wz, with any w; then, the new target particle at z scatters simultaneously with the two 
projectile particles at x and y. 

We see that the two-particle correlation {t^ty) — (t^) {ty) gets built from 'target fluctuations' 
(i.e., splitting processes leading to a change in the particle number in the target distribution), 
via multiple scattering. This mechanism is quite similar to that identified in the context of 
QCD in Ref. [2], and which has inspired an extension of the Balitsky-JIMWLK hierarchy 
known as the 'Pomeron loop' hierarchy [2,20,21]. The general structure of this 'Pomeron 
loop' hierarchy is indeed similar to that for the scattering amplitudes in the toy model (cf. 
Eqs. (I3.21I) - (I3.22I) ). in the sense of including BFKL terms, non-linear terms responsible for 
unitarization, and fluctuation terms which generate correlations in the dilute regime. However, 
some subtle differences persist, which can be partly attributed to real structural differences 
between the dynamics in QCD and that in the toy model, and partly to the different ways how 
perturbation theory is organized in the two cases. It turns out that it is instructive to understand 
these differences in more detail, as they shed new light on the QCD equations themselves. 
This will be discussed at length in the Appendix, where we shall see the two sets of equations 
are essentially equivalent at low density, but they potentially differ from each other in their 
respective extensions towards the high-density regime. In that sense, the lessons drawn from 
the toy model may suggest improvements of the 'Pomeron loop' equations in QCD. 



4 Analytic results: From BFKL growth to diffusive scaling 

In this section we shall present an analytic study of the evolution described by the toy model 
which mimics the corresponding analysis of the high energy evolution in QCD, as standard 
by now in the literature (see especially Refs. [1,2,30,31,33,36]). We shall thus successively 
address the linear approximation (which corresponds to the BFKL equation [7] in QCD), the 
mean field approximation (the toy-model analog of the BK equation [131), and, finally, the 
effects of particle-number fluctuations, that we shall describe by analogs I with the reaction- 
diffusion problem in statistical physics (so like in Ref. [1] for QCD). 



The pertinence of this analogy for the particle model at hand will be confirmed by the numerical 
analysis of this modef in Sect.|5] 
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For the purposes of this analysis, we need to be more specific about the form of the elemen- 
tary particle-particle scattering amplitude t(x\z). Inspired by the analogy with QCD, in which 
this quantity corresponds to the amplitude for dipole-dipole scattering, we shall choose 

t{x\z) = T exp{—\x — z\). (4.1) 

The QCD-origin of this formula can be recognized as follows: after translating to QCD no- 
tations, that is, X = ln(rQ/r^) with r the current dipole size and tq an arbitrary scale of ref- 
erence, Eq. (14.11) becomes equivalent to r(ri|r2) = rr|/r^, where r< = min(ri,r2) and 
r> = max(ri, r2). With the further identification r = the latter is a good approximation to 
the amplitude for the central scattering between two elementary dipoles in QCD. By 'central' 
we mean, as usual, a collision at zero relative impact parameter. It is also understood that the 
amplitude is averaged over the orientations of the two-dimensional vectors Vi and r2. 



4. 1 BFKL evolution in the toy model 



We shall focus on the equations obeyed by the dipole scattering amplitudes and start with 
the linearized version of Eq. (13.211 ). which plays the role of BFKL equation [7] within the toy 
model. With the specific form of t^z given above, this equation reads 

^ = j dzeM-\x~A)tz (4.2) 
(averages are implicitly assumed). Defining the Mellin transform with respect to exp(— x) 

oo 

t(7,r)= / dx exp(7x)t(x,F), (4.3) 



one finds that ^(7, Y) = t{'j,Y =0) exp[x(7) Y] where the 'characteristic function' ^(7) is the 
Mellin transform of the kernel in Eq. (|4.2|) and reads 

Xil) = 7^ + • (4-4) 

1 - 7 1+7 

Then the general solution to Eq. (14.21) is obtained as 

t{x, = j Kl, 0) exp[x(7)l^ - ix] , (4.5) 
c 

where the integration contour C is parallel to the imaginary axis and with its real part being such 
that |Re(7) | < 1. From Eq. (14.31) it is clear that t(7, 0) is the Mellin transform of the scattering 
amplitude at F = 0. For definiteness, we shall assume that the target is initially composed a 
single dipole of size ro = exp(— xo/2), in which case we have ^(7, 0) = rxi'j) exp(7Xo). 

Let us now analyze some special limits of Eq. (14.51) . to which we shall refer by using a 
QCD-inspired terminology: 

(i) The Pomeron intercept : Consider the high-energy limit F ^ 00 at fixed x (which is 
the limit used in QCD to define the 'Pomeron intercept'). In the kinematical plane (F, x), this 
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corresponds to the evolution along a (nearly) vertical axis. In this limit, the integral is dominated 
by the region around the point 7^ that satisfies x'(7p) = 0- One easily finds 7p = and x(7p) = 
2, and by also performing the Gaussian integration around this saddle point one obtains 



r 



{x - Xq) 
8Y 
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t{x,Y) 



exp 2Y - 



(4.6) 



The F-dependence of the amplitude is similar to the QCD one, as in both cases it increases ex- 
ponentially in Y. (Recall that in QCD the dominant F-dependence in this limit is ~ exp(co'pF), 
with Up = 4 In 2.) However the analogy is not complete, since the toy-model amplitude as given 
above is independen{l£] of the projectile dipole size r. By contrast, the respective amplitude in 
QCD is proportional to r, which is however a relatively 'weak' dependence when compared 
to the corresponding result at fixed order perturbation theory, namely t(r) ~ ('color trans- 
parency'). Note also the Gaussian dependence of Eq. (14.61) upon x — xq, which describes diffu- 
sion with a diffusive radius oc F^/^. This is again similar to the well-known BFKL diffusion in 
the logarithmic variable ln(l/r^). 

(ii) Double logarithmic approximation (DLA) : This is the evolution along a direction in the 
plane (Y, x) which is such that the difference x — xq between the dipole sizes increases faster 
than the rapidity Y. When x — Xq ^ Y, the integral is dominated by a value of 7 which is 
close to 1. The saddle point occurs at 7dla = 1 — [Y/ — xq)]^^^ and we are lead to 



This is exactly the same as the corresponding QCD result, which should not come as a surprise: 
in the 'double logarithmic limit' of QCD, one of the daughter dipoles has a size equal to the 
parent one, a feature which is a built-in in the one-dimensional model under consideration. In 
particular, the above result shows 'color transparency', that is, t(x) oc exp[— (x — xq)], in full 
analogy with the QCD behaviour t{r) ~ at very small values of r. 

4.2 Unitarity corrections in the mean field approximation 

As explicit in the previous analysis, the high-energy evolution described by Eq. (14.21 ) leads 
to a scattering amplitude which rises rapidly with Y (cf. Eq. (14.61) ). and thus eventually violates 
the unitarity bound t < 1. So, to study the high-energy behaviour, one needs to include the 
non-linear terms responsible for unitarity corrections. We shall first do so within the mean field 
approximation {txty) ~ {tx) {ty), in which Eq. (13.211) reduces to a closed, non-linear equation 
(the expectation values are again implicit) 



In the QCD terminology, the high-energy solution (14.61 ) exhibits an 'anomalous dimension' 1— 7p = 1, 
which is the maximal possible value; for comparison, this value is 1/2 for the corresponding solution to 
the BFKL equation. 



t{x,Y) 



T 



[(x — xo)F] "^/^exp — (x — Xq) + Y^4(x — Xo)F 



(4.7) 



2v^ 




(4.8) 
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which is the toy model analog of the BK equation in QCD [12, 13]. This analogy is not only 
formal, but it covers several essential aspects: (a) Like the BK equation, Eq. (14.81 ) is consistent 
with unitarity and, moreover, its solution approaches the unitarity bound at large Y (as obvious 
from the fact that t = 1 is a fixed point), (b) In the weak scattering limit t <^ 1, Eq. (14.81) 
reduces, as we have seen, to a linear equation which describes an exponential increase with Y 



and diffusion in x. Together, these properties imply Eq. (14.81) is in the universality class of the 
FKPP equation [4], so like the BK equation itself [33], which in turn implies that the evolution 
towards saturation is driven by the linear dynamics in the dilute regime — the precise form of 
the non-linear terms is unimportant so long as they provide saturation. It is then easy to adapt 
the analysis of the BK equation in Refs. [17, 31, 33] to Eq. (14.81) . with the following results: 

For sufficiently large values of Y (in order to loose memory of the initial condition and 
reach a universal behaviour), the solution t(x, Y) is a traveling wave, that is, a front which 
interpolates between t = 1 at large negative values of x and t at large positive values 
and which propagates towards larger values of x when increasing Y. The position Xs{Y) of 
this front defines the saturation line, that is, the direction of evolution in the plane (Y, x) along 
which the scattering amplitude is constant and of 0{1). As mentioned before, the location of 
this line can be inferred from the solution to the linearized ('BFKL') equation (14.21 ). We thus 
return to this solution, Eq. (14.51 ). and consider a third direction of evolution in the (Y, x)-plane, 
that corresponding to the saturation line: 

(iii) The saturation line : When increasing the rapidity along this line, the position x (i.e., the 
inverse dipole size) should be correspondingly increased in order for the amplitude to remain 
constant. For Y large enough, one can use the saddle point approximation in Eq. (|4.5I) . The 
saddle point condition 

x'{7s)Y -xs = 0, (4.9) 

together with the condition that the exponent vanishes (in order for the amplitude to be roughly 
constant) along the saturation line 

x{ls)Y~^sXs = 0, (4.10) 

uniquely determine the value of the saturation saddle point 7^ and the line XsiY) (the latter up 
to an additive constant), which in this approximation is simply a straight line: Xs{Y) ^ XsY. 
Namely, one finds 

X'(7.) = ^^7. = 4, and A. = ^ = 3v^. (4.11) 
Is V3 75 

It is in fact possible to improve our estimate for As by more properly taking into account the non- 
linear effects. (Within the context of the linear equation, this requires introducing an absorptive 
boundary parallel to the saturation line; see Ref. [31] for details.) One thus finds 

dY -fs 2'ysY 2Y 

in agreement with results from the FKPP equation [4]. To the same accuracy, the amplitude in 
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the vicinity of the saturation line can be obtained via an expansion around the saturation saddle 
point, which yields 



t{x,Y) = cit{x - Xs + C2)exp -'^si.X-Xs) 




(4.13) 



This is strictly valid for 1 ^ a; — ^ 2'x"{^s)Y , which is a parametrically large window at 
large Y . In Eq. (14.131) . Ci and C2 are unknown constants of order 0{\), x"{ls) = 27, and the 
linear factor x — Xg arises from the absorptive boundary. 



We notice that in the region x — Xg <^ J2x"{'~is)Y , where the diffusion term in the exponent 



in Eq. (14.131) can be neglected, the amplitude exhibits geometric scaling, i.e. it depends on x 
and Y only through the combined variable x — Xs{Y). This is what one means by a 'traveling 
wave' [33] : a front which when increasing Y gets simply translated towards larger values of x, 
without being distorted. 

It is furthermore instructive to translate the above picture to the average particle number in 
the target, n{x, Y). In the tail of the front at x XsiY), where t ^ 1, we have t(x, Y) ^ 
/ (izT{x\z)n{z,Y), where the 2;-integration is peaked dX z = x. Thus, a 'small' (large-x) in- 
coming 'dipole' scatters off target 'dipoles' with a similar 'size' {z ~ x), and the weakness 
of the interaction corresponds to the fact that the target looks dilute on this resolution scale: 
n(x, Y) <^ 1/t. When increasing Y in this dilute regime, n rises very fast, exponentially in 
Y, so like t. Around the front position (x ~ XsiY)), one has t ~ 0{l) and thus n ~ 0{1/t): 
on this resolution scale, the target is dense. Finally, behind the front, where t saturates at 1, the 
occupation number does not saturate, rather it keeps growing with Y, albeit only slowly (cf. 
Eq. (13.51 )): n{x,Y) ~ (F — Yc)/t, where Yc ~ ln(l/r) is the 'critical' rapidity for reaching 
saturation in the deposit rate at x : /(x) = 1/r for F > Y^,. The precise value of Yc depends 
upon X and the initial conditions at F = 0. 

4.3 Particle number fluctuations 

As discussed in Sect. 13.21 the evolution equations involve fluctuation terms, like the last 
term in Eq. (13.221) . which reflect the discreteness of the particle number and are inconsistent 
with the mean field approximation underlying Eq. (14.81) . Since the occupation numbers are 
large at saturation, n ~ 1/r ^ 1, one may still hope that the effects of such fluctuations 
are relatively small and can be treated in perturbation theory around the mean-field results. 
However, this expectation turns out to be naive, as only recently understood in the context of 
QCD [1,2, 19,50]. The high-energy evolution is in fact dramatically sensitive to fluctuations, 
and the physical reason for that is implicit in the previous discussion: the evolution is driven by 
the BFKL growth and diffusion in the dilute tail of the front, where the occupation numbers are 
small (of order one) and thus the effects of particle-number fluctuations are important indeed. 

The previous analysis of the linear and the mean-field approximations suggests that the 
model under consideration falls in the universality class of the 'reaction-diffusion' process [4, 
51] (cf . the Introduction), and this conclusion will be further supported by the numerical analysis 
in the next section. It is then possible — like in the corresponding studies in QCD [1,2] — to 
rely on the correspondence with statistical physics in order to characterize the universal aspects 
of the dynamics in the presence of fluctuations. These aspects are mostly qualitative and refer 
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to the 'late-time' (here, large Y) behaviour at 'weak coupling' (here, r ^ 1). In fact, the semi- 
quantitative results which are known to be universal are valid only for values of r which are too 
small to be relevant for the analogy with QCD (recall that r = a^). Yet, most of the qualitative 
features that we shall enumerate below will be later recognized on the numerical results in Sect. 

m 

(i) For a given initial condition at F = 0, the stochastic evolution up to Y generates a sta- 
tistical ensemble of fronts (rather than a single front for the deterministic dynamics in the 
MFA), with the different fronts in the ensemble differing by their respective front positions. 
Hence, XsiY) becomes now a random variable. 
(ii) To a very good approximation, the distribution of Xg at F is a Gaussian, with an expec- 
tation value {xg) and a dispersion cr^(V) = {xD — {xs)"^ which rise both linearly with 
Y: (xg) = XgY and cr^(F) = DfY. The limitations of this Gaussian approximation have 
been analytically studied in Refs. [6,37] in the limit of high occupancy (or weak coupling), 
which amounts to letting r ^ within the toy model, 
(iii) The (asymptotic) value of the average velocity Xg is significantly smaller than its mean- 
field prediction in Eq. (14.1 II) . The respective deviation is analytically known in the formal 
weak-coupling limit [5], with a result which when adapted to the present toy model reads 



A. -^-^4^ = 3^3-^, (4.14) 



2 



valid when In^ r ^ 1. In the same limit, the front diffusion coefficient Df is known to 
behave like Di ~ l/ln^(l/r) [6], and thus to vanish, at it should, when r — > 0. Notice, 
however, the very slow, logarithmic in r, convergence of these results to their respective 
mean-field limits, which reflects the strong sensitivity of the evolution to fluctuations, 
(iv) When increasing Y, the asymptotic value of As is rapidly reached (exponentially in Y), as 
opposed to the rather slow convergence predicted by the mean-field result (14.121) . 
(v) The individual fronts in the ensemble exhibit geometric scaling, but only over a compact 
region: that is, for each front, there is only a finite distance x — Xg < Axg, with Axg ~ 
(l/7s) ln(l/r), ahead of the front where the amplitude scales like t{x) oc e"'''^^^"^'^), with 
7s = 1/ v^. This should be contrasted to the mean-field amplitude in Eq. (|4.13l) . for which 
the scaling window x — Xg (x VY is ever increasing with Y, and hence it can become 
arbitrarily large. 

(vi) The width Axg of the scaling window in the presence of fluctuations is precisely the dis- 
tance over which the amplitude fails from its value t ~ 1 at saturation to a value of 0{t), 
which is the value corresponding to an occupation number n{x) of This is not an ac- 
cident: in an actual event, the occupation number per site is discrete, n{x) = 0, 1, 2, etc., 
so it cannot fall down to a non-zero value which is smaller than one. Hence, a front in n{x) 
has naturally an end point, which is the rightmost occupied bin, and where the occupation 
number is of (9(1). Now, the corresponding front for t(x) is not truly compact, because of 
the non-locality in the elementary scattering amplitude (recall that t(a:;) ^ / dzT(x\z)n{z) 
in the dilute regime). However, the latter decays very fast at large distances, cf. Eq. (|4.1I) . 
so for x — Xg':$> Axg the front t{x) develops an exponential tail of the 'color transparency' 
type: t{x) oc e^^. 

(vii) Physical observables, like the average dipole amplitude (t(x))y, are obtained by averaging 
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over the ensemble. In this averaging, the geometric scaling property of the individual fronts 
is washed out by the dispersion in the positions of the fronts [1, 19], and is eventually 
replaced, at sufficiently large Y, by a new type of scaling [1,2], known as diffusive scaling 
[36]. Specifically, is estimated as 



Y 



-Erfc 
2 



X - (Xs) 



where Erfc(x) is the complimentary error function: 



(4.15) 



Erfc(x) 



2 - 



exp(— x 



exp(— X 



for X ^ — 1 
for |x| ^ 1 
for X ^ 1 . 



(4.16) 



The diffusive scaling approximation, Eq. (14.151) . holds when DfY ^ 1 and in a wide 
region around the (average) saturation line (xs)y, such that |x — (x^) | ^ ^sD{Y. It also 
smoothly interpolates towards the black disk limit (t(x)) = 1 as x — oo. For even larger 
distances x — (x^) ^ 7s-DfF, one recovers the color-transparency behaviour, (t(x)) oc 
exp[— (x — {xs))], with a prefactor which is sensitive to fluctuations (see Refs. [2, 36] for 
details). 



5 Some numerical results 

As mentioned in the Introduction, our present numerical investigation of the one-dimensional 
model is only exploratory, with the purpose of illustrating some of the analytic results and qual- 
itative features anticipated in the previous section and, especially, verifying that the present 
model lies indeed in the universality class of the reaction-diffusion process. In this section, we 
shall briefly present our numerical technique (a Monte-Carlo simulation) and then expose the 
first results obtained in this way. Further verifications and more detailed results are deferred to 
a subsequent publication. 

We would like to simulate the particle evolution described by the master equation (12.41 ) along 
with the deposit rate in Eq. (12.91) . To that aim, we need to specify the one-dimensional lattice 
(that is, the extremal values, Xmin and Xmax, for x, and the lattice spacing A) together with the 
initial conditions at F = and the value of the parameter r which enters the elementary scat- 
tering amplitude (14.11) . In practice, we shall use a lattice of 1 152 size bins between Xmin = — 16 
and Xmax = 128 (hence, A = 0.125). Also, we choose r = 0.01, which via the identification 
T = al corresponds to a value = 0.1 for the coupling constant in QCD. 

We are merely interested in the universal aspects of the evolution at large values of Y, which 
we expect to be largely insensitive to the precise form of the initial conditions. For definiteness, 
we choose the particle occupation number at Y = to be a step function: n(x, Y = 0) = 
Nq9{—x), with Nq = 2. For this distribution, we can compute the initial value of the deposit 
rate fi according to Eqs. (12.91 ) and (14.11 ). 
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Fig. 1. Event-by-event evolution of the scattering amplitude: (a) Rapidity evolution of the scattering 
amplitude for a single event, (b) Amplitude for 10 events (dashed lines) and average amplitude (solid 
line) for Y = 5, 12.5, and 20. 



To evolve the particle system, we have to simulate particle creation according to the deposit 
rate /«. One physical step in this evolution, which consists in the emission of an additional 
particle, involves three steps in the actual algorithm: 

(1) we first randomly select the rapidity SY at which the emission occurs. This is done by 
generating an exponentially decreasing distribution in 'time' (Y) in which the relaxation 
time is fixed as the inverse of the total emission rate J2i /«; 

(2) the site j at which the new particle is created is fixed by randomly selecting a site according 
to a probability law specified by the deposit rate (that is, to each site Xj, we attribute the 
weight function /j); 

(3) the particle number is increased by one at the selected site j and the deposit rates fi are 
updated for all the sites. To speed up this last step, it is useful to notice that, the only 
change in going from fi(Y) to fi{Y + 6Y) refers to the replacement of nj by rij + 1. This 
observation together with Eq. (12.91 ) immediately imply 

MY + 6Y) = -(1 - a,,) + aijUiY). (5.17) 
r 

Those steps are repeated up to a maximal rapidity Y^^y^, that we have chosen as y^ax = 20. To 
study the stochastic aspects of the evolution, we have generated A^ev = 10^ distinct events. 

To begin with, let us consider the rapidity evolution of a single event (see the left plot in 
Fig. [U). This demonstrates the formation of the traveling-wave pattern through two distinct 
mechanisms: 

• particles are produced within already occupied sites. This evolution is basically governed by 
the mean-field dynamics, cf. Eq. (14.81 ), and results in the formation of a traveling front with 
the critical parameters discussed in the previous section; the (average) velocity of this front 
is however influenced by the fluctuations to be discussed below (cf. Eq. (14.141) ). 

• from time to time, a rare fluctuation appears far ahead the tip of the front, in the region that 
was previously unoccupied. (Such fluctuations are rare since their probability is exponen- 
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(a) Number of particle per bin (b) Scattering amplitude 

Fig. 2. Evolution of the average quantities (particle number and scattering amplitude) obtained after 10^ 
events. Results are displayed as a function of x for, from left to right, y = 0, 2, 4, . . . , 20. 



tially decreasing with the distance to the already occupied bins, cf. Eq. (14.11) .) The subsequent 
evolution is then a competition between the local, BFKL-like, growth of that particular fluc- 
tuation and the growth and progression of the mean-field-like wavefront. When considering 
a set of events, such fluctuations lead to an increasing dispersion in the position of the fronts 
(see right plot of Fig. [T]), to which we shall shortly return. 

By averaging over a huge number of events {N^^ = 10^), we obtain the averaged distributions 
shown in Fig. [21 for the average number of particles per lattice site and for the average scattering 
amplitude. In the left plot, we observe large fluctuations in the average number of particle per 
site when (n) ^ 10"^ = l/Ney. Those correspond to rare fluctuations (only a few events 
have nonzero occupation number in the respective bins) in the dilute tail of the front. It is also 
interesting to notice that the fluctuations start to be visible around {n) = 1/ \/N^ f» 0.003, as 
expected for particle-number fluctuations. One can also check that, within the dense regime, 
the number of particles per bin increases linearly as expected from Eq. (13.51) . 

Concerning the average amplitude, the fact that this remains smooth even in the dilute tail of 
the front is just because this quantity is obtained by a convolution of {n) with the elementary 
interaction Tij (recall that (tj) Y.j 'Tij {nj) in the dilute regime). Still for the average amplitude, 
although one can notice a traveling-wave pattern in the right plot in Fig. [21 a closer look reveals 
that the shape of the average front is in fact changing with Y — its slope decreases with Y. This 
reflects the violation of geometric scaling through fluctuations, an important physical effect that 
we shall now study in more detail. 

To that aim, we need to first study the statistics of the position of the front ('the saturation 
scale Xs') as a function of rapidity. For a given event t(x, Y), the position of the front is deter- 
mined by t{x = Xs{Y), Y) = tQ, with to some constant number that we have here chosen as 
to = 0.1 (but we checked that different choices do not significantly modify the results). From 
the discussion in Sect. [H we expect both the average value of Xs{Y) and the associated disper- 
sion to grow linearly with Y. For the average value, this is checked on the left plot in Fig. [3l 
which exhibits our numerical results for dy {xs{Y)). This plot further shows that the large-F 
asymptotic of the velocity is attained pretty fast and, moreover, this asymptotic value is sig- 
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(a) The average velocity of the travehng front 




(b) Dispersion in the position of the fronts 



Fig. 3. The statistics of the front position (or 'saturation scale') Xg. Both the average position and the 
squared dispersion increase linearly with rapidity. 




-2 2 

{x-Xs{Y))/a{Y) 

Fig. 4. The average amplitude represented as a function of the diffusive scaling variable. When rapidity 
increases, the asymptotic result (14.151) (indicated in the figure by 'limit') is indeed obtained. 



nificantly smaller than the corresponding prediction of the mean-field approximation, namely 
Xs = 3v^ (cf. Eq. (|4.1 II) ). Next, we turn the corresponding dispersion, cr^ = (s^) — (s^)^, with 
the results displayed in the right plot in Fig.[3l This is well fitted by a linear increase, ~ DiY , 
with a 'front diffusion' coefficient = 1.224. Thus, the statistics of the 'saturation scale' fol- 
lows indeed the pattern expected for a 'reaction-diffusion' process. We are not in a position 
here to also check more quantitative aspects, so like the slow approach towards the respective 
mean-field results with decreasing r, since the parameter r in our calculation not only has a 
fixed value, but this value is also too large for analytic formulae like Eq. (|4.14l) to apply. 

We now return to an analysis of the shape of the average amplitude (cf. the right plot in 
Fig. |2l), with the purpose of verifying another important prediction of the 'reaction-diffusion' 
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problem: the emergence of diffusive scaling at large Y (cf. Eq. (14.151) ). To check that, we have 
plotted the average amplitude obtained from our numerical simulations as a function of the 
expected scaling variable (x — {xs))/a, with (xg) and taken from our previous analysis (Fig. 
[3]) at various values of Y. As manifest on Fig. |4l the measured average amplitude converges 
indeed towards the expected asymptotic curve, Eq. (I4.15|) . when increasing Y. The diffusive 
scaling property is thus fully satisfied by our model. 

6 Conclusions 

In this paper we have presented a (l+l)-dimensional stochastic particle model which mimics 
high-energy evolution and scattering in QCD at fixed impact parameter. The 'time' coordinate 
Y in this model corresponds to rapidity, while the 'spatial' coordinate x represents the (loga- 
rithm of the) transverse size of a dipole in QCD. The model is interesting in several respects: 

(i) The structure of the model is consistent with, and to a large extend determined by, general 
physical principles which are known to hold in QCD: boost-invariance, multiple scatter- 
ing, and evolution via the emission of an additional particle per unit rapidity. 

(ii) Related to the above, the model exhibits a saturation mechanism which is similar to gluon 
saturation in QCD: The new particles created by the evolution are coherently emitted from 
the preexisting ones, with an emission rate which saturates at high density because of the 
multiple scattering between the additional particle and its 'parents', 
(iii) The model appears to be in the universality class of the reaction-diffusion process, as also 
expected for QCD, but unlike other related models in the literature, it does not involve 
explicit vertices for particle recombination. These would be redundant since saturation is 
anyway ensured by coherence effects in the particle emission, as explained above. 

(iv) The model exhibits all the qualitative features expected in QCD at fixed impact param- 
eter, concerning both the mean-field aspects (a BFKL-like growth in the dilute regime, 
the formation of a saturation front with geometric scaling) and the effects of fluctuations 
(dispersion in the positions of the fronts, breakdown of geometric scaling in the statistical 
description, convergence towards diffusive scaling at high energy), 
(v) The model is simple enough to allow for detailed numerical investigations. Hence, it can 
be used as a playground to develop and test our intuition concerning the dynamics in QCD 
at high energy. Complex phenomena which are expected, but not yet demonstrated, in 
QCD (like the emergence of diffusive scaling) can be visualized in this context and their 
physical consequences can be explicitly analyzed. 

(vi) The structural aspects of the model are particularly interesting since they may inspire our 
searches for corresponding structures in QCD. In particular, the evolution equations for the 
scattering amplitudes in this model appear as a natural generalization (within the limits of 
the model, of course) of the Balitsky-JIMWLK equations in which the projectile and target 
are symmetrically treated: for each of these two systems, the equations include particle- 
number fluctuations and multiple scattering with both the particles internal to the system 
(saturation effects) and those in the other system (unitarity corrections). 

In the remaining part of this concluding section we would like to elaborate a little bit more on 
the last point above and comparatively discuss the evolution equations in the toy model and the 
'Pomeron loop' equations proposed in QCD at large [2,20,21]. This comparison is discussed 
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in more detail in the Appendix from which we shall extract here the main conclusions. 

By inspection of the two sets of equations — see, e.g., Eqs. (13.211 ), (13.221 ) and (lA.ll) for 
the toy model, and, respectively, Eqs. (2.7)-(2.8) in Ref. [21] for QCD — , one can notice 
some important differences which formally refer to the structure and physical interpretation 
of the respective fluctuation terms: Within the toy model, we have seen (cf. Eq. (13.231) ) that 
these terms correspond to the multiple scattering of only one among the two child particles 
produced by a splitting in the target with two (or more) particles from the projectile. By contrast, 
these particular multiple-scattering effects are not at all included in the respective 'Pomeron 
loop' equations, where the fluctuation terms rather refer to the separate scattering of both child 
dipoles produced after a splitting in the target with two (different) dipoles in the projectile. 
However, as explained in the Appendix, such differences seem to be inessential (except in the 
very early stages of the evolution) since they amount to a reorganization of the perturbation 
theory in the dilute regime, which however leads to the same dominant behaviour at sufficiently 
high energy — namely, such that F ^ 1 in the toy model and, respectively, cisY ^ I'm QCD. 

To briefly explain this reorganization, let us consider the weak-scattering regime, where the 
scattering amplitude for a projectile made with k particles can be estimated as (for F ^ 1) 



lt^^\xi,...,Xk))^ = J T{xi\zi)---T{xk\zk){n'-^\zi,...,Zk))^, (6.1) 



where {n^'^^) is the 'normal-ordered' /c-body density of Eq. (13.61) . This formula assumes that 
the k external particles scatter off k dijferent particles in the target, and thus treats the target and 
the projectile on the same footing — multiple scattering is neglected for both of them. Strictly 
speaking, Eq. (|6.1I) is not the same as the fc-particle scattering amplitude in fixed order per- 
turbation theory (the latter would include additional contributions associated with the multiple 
scattering of individual target particles; see, e.g., Eq. (IA.3I) in the Appendix), but it becomes 
equivalent to the latter when Y ^ 1 since it captures the contribution which has the fastest rise 
with Y (the A;-Pomerons piece of the complete amplitude). 

In QCD at low density (and large Nc), one can define 'bare' amplitudes similar to Eq. (16. lb 
in terms of dipoles and then deduce linear evolution equations for these quantities within the 
dipole picture [21]. The 'Pomeron loop' equations are then obtained by completing these linear 
equations with the non-linear terms responsible for unitarity corrections, as taken over from the 
Balitsky-JIMWLK equations. This procedure amounts to dressing the 'bare' amplitudes {t^^'^) 
with multiple-scattering effects on the side of the projectile. In the toy model, on the other 
hand, the evolution equations obtained in Sect. 13.21 correspond to dressing {t^'^^) with multiple 
scattering on both the projectile, and the target, sides. 

Clearly, the last procedure is conceptually more satisfactory, as it provides a symmetric treat- 
ment of the target and the projectile in the general, strong-scattering, regime. In QCD, such a 
symmetric description is still lacking and its construction remains as an important open prob- 
lem. But although conceptually, and also aesthetically, more appealing, it is not clear to us 
whether such a symmetric description would differ indeed — in terms of physical consequences 
at high energy — from that provided by the already known 'Pomeron loop' equations. 
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A On the fluctuation terms in the Pomeron loop equations 

In this Appendix we shall describe and clarify some structural differences between the evo- 
lution equations generated by the toy model and the Pomeron loop equations in QCD [2,20,21]. 
To that aim, it is useful to also have under ones eyes the third equation in the toy-model hierar- 
chy, as obeyed by {txtytz); this reads: 



d \txtytz) f T~xw /, , , /-I J. \\ , 

— — = J — K^wtytzO- -^x)) + permutations 

to 

j xw wy |/^^^^^]^ _ _ i^)^ _)_ permutations 



+ / (t4i-t,)(i-g(i-t,)). (A.i) 



Already a superficial comparison between the toy-model equations, cf. Eqs. (I3.21l) - (l3.22h 
and (lA.ll) . and the corresponding equations in QCD, cf. Ref. [21], reveals some important dis- 
similarities that we summarize here: 

(A) The equation for {t^) = {tx^ ■ ■ - tx^) within the toy model involve several types of fluc- 
tuation, with the following generic structures: r(t'^~^), r^(t'^^^), . . . , r^^^(t). These terms are 
all of the same order in the low density regime where {t) ~ r, and hence contribute on equal 
footing to building up the many-body correlations. By contrast, the corresponding!^ 'Pomeron 
loop' equation for (T^*^)) involves just one fluctuation term, with the generic structure t{T^^^^'>) 
(where r ~ 0{al) in QCD). 

(B) As discussed in relation with Eq. (13.231) . the fluctuation term in the equation for {txty) 
corresponds to the double scattering of one of the daughter particles produced by a splitting in 
the target. A similar conclusion applies to the higher equations in the toy-model hierarchy, like 
Eq. (lA.ll) : a term like (rxw^wy/T) {tw^z) in the second line of Eq. (lA.ll) describes a process in 
which the two external particles at x and y scatter off the daughter particle produced at w after 
a splitting in the target (with the third external particles at z being just a spectator); also, a term 
like {TxwTywTzw / t) {tw) in the third line there corresponds to a triple scattering for the target 
particle at w, which now scatters off all the external particles. 

By contrast, in the construction [2, 20, 21] of the fluctuation terms in the Pomeron loop 
equations in QCD, one has completely neglected the multiple scattering for the individual target 
dipoles. Rather, the fluctuation terms included in these equations describe the separate scattering 



Note that we use the upper-case notation T for the scattering amplitude of a dipole in QCD to distin- 
guish from the respective amplitudes t in the toy model. 
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of both child dipoles produced after a splitting in the target with a pair of projectile dipoles. 

The above discussion rises several questions: 

(i) Why there are no fluctuation terms in the toy-model equations corresponding to the single 
scattering of both particles w and z produced after a splitting w wz in the target ? 
(ii) Was it correct to neglect multiple scattering for the individual target dipoles in the corre- 
sponding equations in QCD ? 
(iii) What should be the complete structure of the fluctuation terms in QCD ? 

We are not so ambitious to try and fully answer question (iii) in what follows, but we 
would like to clarify at least the answers to the first two questions. Specifically, we will show 
that the reason why the two sets of equations look so different is because they apply to different 
quantities which however have the same dominant behaviour at high energy (namely, for Y ^ 
1). 

The following considerations will be restricted to a dilute target, off which the individual 
projectile particles scatter only once. It is then preferable to express the scattering amplitudes 
in terms of particle densities in the target and follow the evolution of the latter. By using 



t{x) = 1 — exp 



dz n(z) In a(x\z) 



dzr{x\z)n{z) (A. 2) 



(with the last, approximate, equality valid in the dilute regime), together with the definition 
(13.101 ) of the normal-ordered pair density, one can immediately deduce 

(tx^y) ~ j '^xz'Tyw {l^z'^w) J '^xz'Tyw ^'^iji^ ~l~ j '^xz'^yz iP'z) i (A. 3) 

z,w z,w z 

where the contribution proportional to (n*^^)) describes the scattering between the two projec- 
tile particles and two different target particles, while that involving (n) describes the double 
scattering of a same target particle. We shall simultaneously consider the following quantity 

i'xt) = I rxzryu,U2), (A.4) 



xy I ~ j 'xz'yw yi'zw I 1 
z,w 

which represents the toy-model analog of the 2-dipole scattering amplitude (T^^^) which enters 
the 'Pomeron loop' equations in QCD [2,21]. As compared to Eq. (IA.3I) . the double-scattering 
terms are neglected in Eq. (IA.4I) . which therefore treats more symmetrically the target and the 
projectile in this weak-scattering regime (in the sense that all the particles undergo single scat- 
tering). However, it is a priori the less symmetric formula, Eq. (|A.3I) . which yields the 2-particle 
scattering amplitude at fixed order in perturbation theory: indeed, both terms in the r.h.s. of 
Eq. (IA.3I) are of C(t^). A similar discussion applies to the A;-particle amplitudes — (t'') and 
respectively {t^'^^) — with A; > 3. 

In what follow we shall use the evolution equations obeyed by the particle densities in the 
dilute regime to deduce the respective equations for (t^) and (t'-^-'). For the former, we shall 
recover, as expected, the linear part of the general equation (13.221) . but with a more transparent 
interpretation for the fluctuation terms, which will also answer question (i) above. For the latter, 
we shall find the toy-model analog of the second equation in the 'Pomeron loop' hierarchy. 
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The low-density versions of Eqs. (13.121 ) and (13.131 ) read 



(A.5) 



dY 



u 




(A.6) 



dY 



T 



u 



Notice the terms linear in (n) in the r.h.s. of Eq. (IA.6I) : these are fluctuation terms which generate 
a pair (2;, w) of particles in the target via the splitting of a single original particle at z, or at w. 
These terms are very similar to those generated by the dipole picture in QCD [2]. 

The equation obeyed by (txty ) can now be obtained by taking a derivative w.r.t. Y in Eq. (IA.3I) 
and then using Eqs. (IA.5I) - (IA.6I) . By inspection of these equations, one may expect this proce- 
dure to produce two types of 'fluctuation terms' (i.e., terms linear in (t)) : those generated by 
the BFKL evolution of the average particle density, cf. Eq. (|A.5I) . and the 'genuine' fluctuation 
terms coming from the evolution of the pair density, cf. Eq. (IA.6I) . The latter would correspond 
to the separate scattering of both particles w and z produced by the splitting in the target, and 
thus would be the toy-model analog of the fluctuation terms which appear in the Pomeron loops 
equations in QCD. However, from the previous discussion of Eq. (13.231) . we know already that 
the relevant fluctuation term comes fully from the BFKL evolution of (n), cf. Eq. (IA.5I) ! What 
happens then to the fluctuation terms which are explicit in Eq. (IA.6I ) ? As we explain now, they 
are in fact absorbed — via Eq. (IA.3I) — in the structure of the 2-particle amplitudes which enter 
the 'BFKL' terms in Eq. (13.221) . Indeed, Eqs. (IA.3I) and (IA.6I) can be used to deduce 



where in writing the first equality we have already reabsorbed the fluctuation terms into 2-body 
densities like {riunj}, via Eq. (13.101) . In the second line of the above equation one can recognize 
the BFKL terms from Eq. (|3.22l) . as anticipated. 

This answers question (i) above — the analogs of the fluctuation terms in the 'Pomeron 
loop' equations are implicitly included in the toy-model equations, as multiple-scattering con- 
tributions to the amplitudes which enter the BFKL terms there — but rises new questions in 
return : Why are these terms explicit in the QCD equations in Ref. [2, 20, 21], while they are 
implicit in the corresponding toy-model equations ? And why the fluctuations terms which are 
explicit in the toy model (those corresponding to multiple scattering off a same target particle) 
appear to be absent in QCD ? 

The answer to these last questions, and also to the original question (ii), is that the 'Pomeron 
loop' equations are in fact written for different quantities, namely the 'reduced' scattering am- 
plitudes {t^''^) (see Eq. (IA.4I) ) which in this dilute regime neglect multiple scattering for both 
the target particles and the projectile ones. This amounts to a different reorganization of the per- 
turbation theory in the weak-scattering regime, which however leads to the same high-energy 
behaviour, as previously argued in Refs. [52,53]. (The subsequent discussion is in fact adapted 




(A.7) 



2 
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from Ref. [53] .) To see this, let us consider the relative importance of the two terms in Eq. (IA.3I ). 
For more clarity, assume that the target starts at F = as a single particle located at xq. Hence 
(n^)Q = 6{z — xo) and (nf2)^ = 0- Then, clearly, the term proportional to (n) in Eq. (|A.3I) will 
dominate in the early stages of the evolution (F <^ 1), until a non-vanishing pair density gets 
first created via fluctuations, cf. Eq. (IA.6I) . But once the latter becomes non-zero, its subsequent 
evolution is much faster — it rises as a 'double-Pomeron', {n^"^^) ~ exp(2c<jpF) with cuf = 2, 
while the single-density rises only as (n) ~ exp(c<jpF) — and thus it dominates the 2-particle 
amplitude in Eq. (IA.3I) for all rapidities Y > 1. 

To summarize, in the high-energy regime where Y is relatively large, Y > 1, but not too 
large (so that we are still in the weak scattering regime; this condition requires Y <^ ln(l/r)), 
the 2-particle scattering amplitude can be estimated simply as (t^^-*), cf. Eq. (IA.4I) . Equivalently, 
this can be obtained by solving the evolution equation for {t'^'^^) (see below) with the initial con- 
dition (t'^^^)o = 0. Note that this initial condition is different from that to be used in relation with 
the original equation (13.221) . which follows from Eq. (IA.3I) and reads {txty)^ = r(x|xo)'r(y|xo). 

As quite clear by inspection of Eqs. (IA.4I) and (|A.6I) . the fluctuation terms in the equation 
for (t^^^ ^ come fully from the process in which the two external particles scatter off both child 
particles produced by a splitting in the target. Specifically, one finds 



dY 



^(2) 



Ml /+(2) 
7- 



dY 



(A.8) 



fluct 



with the following expression for the fluctuation term: 



dY 



fluct 



z,w 

j {^xz'^yw 



'Txw 'Tyz 



i^z) 



' WZ —I 



(A.9) 



where in the second line we have expressed the average particle density in the target (n^) in 
terms of the one-particle scattering amplitude (t„) by inverting the relation {ty) = Tuz {^z), 
valid in the dilute regime. The inverse r^y of the elementary scattering amplitude exists indeed, 
since r^y is positive-definite; in particular, for the function r^y in Eq. (14.11 ). one has 



r 



\x\y) = ^{l-dl)6{x-y). 



(A. 10) 



As anticipated, Eq. (IA.9I) is the toy-model analog of the fluctuation term in the corresponding 
'Pomeron loop' equation (see Eq. (2.7) in Ref. [21]). This analogy extends, of course, to the 
higher equations in the respective hierarchies. 

We have previously noticed that, for F ^ 1, there is no important loss of accuracy if the so- 
lution (t^^^'^ to Eqs. (|A.8I) - (|A.9I ) with initial condition {t^'^'^)o = is used instead of the solution 
(txty) to the more general equation (|3.22l) with initial condition {txty)^ = T(x\xo)T(y\xo). By 
the same token we deduce that, within QCD, the solutions to the Pomeron loop equations with 
initial conditions {T^''^)o = for A; > 2 should provide the correct physical behaviour at high 
energy. This conclusion is further supported by the arguments in Ref. [52, 53]. 
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